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Abstract 


The  Wigner  distribution  is  a  signal  transformation  which  has  its  origin  in  quantum 
mechanics.  It  possesses  some  important  properties  which  make  it  very  attractive  for 
time-frequency  signal  analysis.  The  Wigner  distribution  was  originally  defined  for 
continuous-time  signals.  A  discrete-time  version  of  it  was  proposed  recently. 
Unfortunately,  this  discrete-time  Wigner  distribution  suffers  from  aliasing  effects, 
which  prevent  several  of  the  properties  of  the  continuous-time  Wigner  distribution  from 
carrying  over  straightforwardly.  In  this  report,  a  discrete-time  Wigner  distribution 
which  does  not  suffer  from  aliasing  is  introduced.  It  is  essentially  an  augmentation  of 
the  previous  version,  incorporating  new  information  about  the  signal  not  contained  in 
the  previous  version.  A  variety  of  methods  have  been  proposed  for  the  detection  of  a 
signal,  with  unknown  signal  parameters,  in  a  noisy  envimoment.  In  this  report,  the 
noise  statistics  are  incorporated  to  reveal  that  certain  processing  of  the  Wigner 
distribution  signal  representation  can  lead  to  an  optimal,  and  often  easy  to  compute, 
detection  scheme.  For  the  special  case  of  discrete-time  chirp  signals  in  complex  white 
Gaussian  noise,  it  is  shown  that  the  optimal  detector  is  equivalent  to  integrating  the 
Wigner  distribution  along  the  line  of  instantaneous  frequency.  If  the  position  and  sweep 
rate  of  the  linear  chirp  are  unknown,  then  a  Generalized  Likelihood  Ratio  Test  (GLRT) 
leads  one  to  integrate  the  Wigner  distribution  along  all  possible  lines  in  the 
time-frequency  plane  and  choose  the  largest  integrated  value  for  comparison  to  a 
threshold.  Simulation  examples  of  the  Wigner  distribution  scheme  are  given  to 
demonstrate  the  utility  concerning  the  detection  of  the  proposed  method. 
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Resume 


La  distribution  de  Wigner  est  un  outil  de  mise  en  forme  de  signal  qui  decoule  de  la 
mecanique  quantique.  Elle  possede  certaines  proprietes  importantes  qui  en  font  un  outil 
tres  interessant  pour  l’analyse  temps-frequence  des  signaux.  La  distribution  de  Wigner 
a  ete  definie  a  l’origine  pour  les  signaux  a  temps  continu.  Une  version  a  temps  discret  a 
ete  proposee  recemment.  Malheureusement,  cette  derniere  presente  des  effets  de 
repliement  de  spectre,  de  sorte  que  plusieurs  des  proprietes  de  la  distribution  de  Wigner 
a  temps  continu  ne  peuvent  lui  etre  transferees  directement.  Dans  le  present  rapport,  on 
traite  d’une  distribution  de  Wigner  a  temps  discret  qui  ne  comporte  pas  de  repliement 
de  spectre.  Cette  distribution  consiste  essentiellement  en  une  extension  de  la  version 
precedente,  a  laquelle  ont  ete  ajoutees  de  nouvelles  donnees  relatives  au  signal. 

Diverses  methodes  ont  ete  proposees  pour  la  detection  d’un  signal,  dont  des  parametres 
sont  inconnus,  dans  un  environnement  bruite.  Dans  le  present  rapport,  les  statistiques 
sur  le  bruit  sont  incorporees  afin  de  montrer  que  certains  traitements  de  la 
representation  du  signal  sous  forme  d’une  distribution  de  Wigner  peuvent  mener  a  un 
plan  de  detection  optimal,  comportant  souvent  des  calculs  relativement  simples.  Dans 
le  cas  particulier  des  impulsions  modulees  en  frequence  a  temps  discret  en  presence  de 
bruit  blanc  gaussien  complexe,  on  montre  que  le  detecteur  optimal  est  obtenu  en 
integrant  la  distribution  de  Wigner  suivant  la  ligne  de  frequence  instantanee.  Si  la 
position  et  la  vitesse  de  balayage  du  signal  a  modulation  de  frequence  lineaire  sont 
inconnues,  un  test  generalise  du  rapport  des  vraisemblances  (TGRV)  permet 
Integration  de  la  distribution  de  Wigner  suivant  toutes  les  lignes  possibles  du  plan 
temps-frequence  et  la  selection  de  la  plus  grande  valeur  integree  pour  fins  de 
comparaison  a  un  seuil.  Des  exemples  de  simulation  de  plan  de  distribution  de  Wigner 
sont  donnes  dans  le  but  de  demontrer  l’utilite  en  ce  qui  concerne  la  detection  de  la 
methode  proposee. 
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Executive  summary 


One  of  the  central  problems  in  High  Frequency  (HF)  radar  data  is  the  analysis  of  a  time 
series.  The  Fourier  transform  method,  or  Doppler  processing  method,  has  been 
generally  used  in  HF  radar  to  detect  targets  that  are  moving  with  constant  radial 
acceleration.  Examples  of  accelerating  targets  are  manoeuvring  aircrafts  and  missiles. 
Our  previous  studies  show  that  there  are  limitations  and  shortcomings  in  the  Fourier 
transform  method  to  detect  accelerating  targets  because  of  the  phenomenon  known  as 
Doppler  smearing  [1,2].  Our  studies  also  show  that  when  the  target  is  constantly 
accelerating,  the  Fourier  method  may  still  be  used  to  detect  the  target  and  estimate  its 
median  velocity,  provided  the  acceleration  is  small  enough  in  the  sense  to  be  described 
in  this  report.  It  is  shown  that  for  a  given  acceleration,  the  number  of  pulses  cannot  be 
increased  indefinitely  without  resulting  in  catastrophic  failure  of  the  method. 
Conversely,  for  a  given  number  of  pulses,  the  acceleration  cannot  be  arbitrarily  large 
without  resulting  in  catastrophic  failure  of  the  method.  Thus  the  number  of  pulses  and 
the  acceleration  have  to  be  matched  to  achieve  optimum  performance  [1,2]. 

Consequently,  for  the  interpretation  of  radar  data  in  terms  of  a  changing  frequency 
content,  we  need  a  representation  of  our  data  as  a  function  of  both  time  and  frequency. 
The  purpose  of  this  report  is  therefore  to  stress  the  importance  of  alternative  methods 
which  have  been  given  little  attention  in  the  past,  namely,  the  joint  time- frequency 
representation  of  signals.  In  this  report  we  choose  the  Wigner  distribution  together  with 
a  Radon  transform  as  the  basic  tool,  since  both  possess  some  important  properties 
which  are  very  attractive  for  time-frequency  signal  analysis  [3]. 

A  variety  of  methods  have  been  proposed  for  the  detection  of  a  signal,  with  unknown 
signal  parameters,  in  a  noisy  envirnoment.  In  this  report,  the  noise  statistics  are 
incorporated  to  reveal  that  certain  processing  of  the  Wigner  distribution  signal 
representation  can  lead  to  an  optimal,  and  often  easy  to  compute,  detection  scheme.  For 
the  special  case  of  discrete-time  chirp  signals  in  complex  white  Gaussian  noise,  it  is 
shown  that  the  optimal  detector  is  equivalent  to  integrating  the  Wigner  distribution 
along  the  line  of  instantaneous  frequency.  If  the  position  and  sweep  rate  of  the  linear 
chirp  are  unknown,  then  a  Generalized  Likelihood  Ratio  Test  (GLRT)  leads  one  to 
integrate  the  Wigner  distribution  along  all  possible  lines  in  the  time-frequency  plane 
and  choose  the  largest  integrated  value  for  comparison  to  a  threshold.  Simulation 
examples  of  the  Wigner  distribution  scheme  are  given  to  demonstrate  the  utility 
concerning  the  detection  of  the  proposed  method. 

In  this  report,  a  new  discrete-time  Wigner  distribution  has  also  been  proposed.  The 
ideas  presented  can  be  extented  to  formulate  a  Wigner  distribution  in  both  time  and 
frequency,  for  discrete-time  signals  of  finite  duration,  in  the  same  spirit  as  development 
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of  the  discrete  Fourier  transform.  Application  of  these  ideas  to  actual  data  analysis  is 
the  subject  of  further  research. 


T.  Thayaparan,  A.  Yasotharan.  2001 .  A  Novel  Approach  for  the  Wigner  Distribution 
Formulation  of  the  Optimum  Detection  Problem  for  a  Discrete-Time  Chirp  Signal. 
DREO  TM  2001-141.  Defence  Research  Establishment  Ottawa. 
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Sommaire 


Un  des  principaux  problemes  eprouve  avec  les  donnees  de  radar  haute  frequence  (HF) 
est  1’ analyse  d’une  serie  chronologique.  La  methode  a  transformee  de  Fourier,  ou  la 
methode  de  traitement  Doppler,  a  generalement  ete  utilisee  dans  le  radar  HF  pour  la 
detection  des  cibles  a  acceleration  radiale  constante.  Les  cibles  presentant  une 
acceleration  sont,  par  exemple,  les  aeronefs  qui  effectuent  des  manuvres  et  les  missiles. 
Nos  eetudes  preeceedantes  montre  que  la  methode  a  transformee  de  Fourier  comporte 
des  limites  et  des  lacunes  pour  la  detection  des  cibles  a  acceleration,  en  raison  du 
phenomene  appele  la  persistance  Doppler  [1,2].  Nos  eetudes  preeceedantes  montre 
aussi  que  dans  le  cas  d’une  cible  e  acceleration  constante,  la  methode  de  Fourier  peut 
encore  etre  appliquee  pour  detecter  la  cible  et  estimer  sa  vitesse  mediane,  a  condition 
que  V acceleration  soit  assez  faible  dans  le  contexte  du  present  rapport.  On  montre  que 
pour  une  acceleration  donnee,  le  nombre  d’ impulsions  ne  peut  pas  etre  augmente 
indefiniment  sans  qu’il  en  resulte  un  echec  catastrophique  de  la  methode.  Inversement, 
pour  un  nombre  donne  d’impulsions,  V  acceleration  ne  peut  pas  augmenter 
arbitrairement  sans  qu’il  en  resulte  aussi  un  echec  catastrophique  de  la  methode.  Ainsi, 
le  nombre  d’impulsions  et  1’ acceleration  doivent  etre  adaptes  pour  permettre  des 
performances  optimales  [1,2]. 

Par  consequent,  pour  interpreter  les  donnees  radar  du  point  de  vue  d’un  contenu  en 
frequences  variable,  nous  devons  les  representer  aussi  bien  en  fonction  du  temps  qu’en 
fonction  de  la  frequence.  L’objet  du  present  rapport  est  done  de  souligner  l’importance 
des  methodes  de  rechange  auxquelles  on  a  accorde  peu  d’attention  par  le  passe, 
notamment  la  representation  temps-frequence  combinee  des  signaux.  Dans  le  rapport, 
nous  choisissons  la  distribution  de  Wigner  avec  une  transformee  de  Radon  comme  outil 
de  base,  etant  donne  que  ces  deux  elements  offrent  certaines  proprietes  importantes  qui 
presentent  un  grand  interet  pour  1’ analyse  temps-frequence  des  signaux  [3]. 

Diverses  methodes  ont  ete  proposees  pour  la  detection  d’un  signal,  dont  des  parametres 
sont  inconnus,  dans  un  environnement  bruite.  Dans  le  present  rapport,  les  statistiques 
sur  le  bruit  sont  incorporees  afin  de  montrer  que  certains  traitements  de  la 
representation  du  signal  sous  forme  d’une  distribution  de  Wigner  peuvent  mener  a  un 
plan  de  detection  optimal,  comportant  souvent  des  calculs  relativement  simples.  Dans 
le  cas  particulier  des  impulsions  modulees  en  frequence  a  temps  discret  en  presence  de 
bruit  blanc  gaussien  complexe,  on  montre  que  le  detecteur  optimal  est  obtenu  en 
integrant  la  distribution  de  Wigner  suivant  la  ligne  de  frequence  instantanee.  Si  la 
position  et  la  vitesse  de  balayage  du  signal  a  modulation  de  frequence  lineaire  sont 
inconnues,  un  test  generalise  du  rapport  des  vraisemblances  (TGRV)  permet 
F  integration  de  la  distribution  de  Wigner  suivant  toutes  les  lignes  possibles  du  plan 
temps-frequence  et  la  selection  de  la  plus  grande  valeur  integree  pour  fins  de 
comparaison  a  un  seuil.  Des  exemples  de  simulation  de  plan  de  distribution  de  Wigner 
sont  donnes  dans  le  but  de  demontrer  l’utilite  en  ce  qui  conceme  la  detection  de  la 
methode  proposee. 
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Dans  le  present  rapport,  une  nouvelle  distribution  de  Wigner  a  temps  discret  a  aussi  ete 
proposee.  Les  notions  presentees  peuvent  etre  etendues  a  la  formulation  d’une 
distribution  de  Wigner  tant  dans  le  domaine  temporel  que  dans  le  domaine  frequentiel, 
pour  des  signaux  a  temps  discret  de  duree  finie,  dans  le  meme  esprit  qui  a  ete  adopte 
pour  Elaboration  de  la  transformee  de  Fourier  discrete.  L’ application  de  ces  notions  a 
E  analyse  de  donnees  reelles  fait  Fobjet  de  recherches  plus  approfondies. 

T.  Thayaparan,  A.  Yasotharan.  2001.  Nouvelle  approche  pour  la  formulation  de  la  distribution 
de  Wigner  du  probleme  de  detection  optimale  d’un  signal  a  temps  discret. 

DREO  TM  2001-141.  Centre  de  Recherches  pour  la  Defense  Ottawa. 
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1.  Introduction 


The  Fourier  Transform  is  at  the  heart  of  a  wide  range  of  techniques  that  are  used  in 
radar  data  analysis  and  processing.  Mapping  the  data  into  the  temporal  frequency 
domain  is  an  effective  way  of  recording  the  data  such  that  their  global  characteristics 
can  be  assessed.  However,  the  change  of  frequency  content  with  time  is  one  of  the  main 
features  we  observe  in  HF  radar  data.  Because  of  this  change  of  frequency  content  with 
time,  radar  signals  belong  to  the  class  of  non-stationary  signals. 

One  of  the  central  problems  in  High  Frequency  (HF)  radar  data  is  the  analysis  of  a  time 
series.  The  Fourier  transform  method,  or  Doppler  processing  method,  has  been 
generally  used  in  HF  radar  to  detect  targets  that  are  moving  with  constant  radial 
acceleration.  Examples  of  accelerating  targets  are  manoeuvring  aircrafts  and  missiles. 
Our  previous  studies  show  that  there  are  limitations  and  shortcomings  in  the  Fourier 
transform  method  to  detect  accelerating  targets  because  of  the  phenomenon  known  as 
Doppler  smearing  [1,2].  Our  studies  also  show  that  when  the  target  is  constantly 
accelerating,  the  Fourier  method  may  still  be  used  to  detect  a  target  and  estimate  its 
median  velocity,  provided  the  acceleration  is  small  enough  in  the  sense  to  be  described 
in  this  report.  It  is  shown  that  for  a  given  acceleration,  the  number  of  pulses  cannot  be 
increased  indefinitely  without  resulting  in  catastrophic  failure  of  the  method. 
Conversely,  for  a  given  number  of  pulses,  the  acceleration  cannot  be  arbitrarily  large 
without  resulting  in  catastrophic  failure  of  the  method.  Thus  the  number  of  pulses  and 
the  acceleration  have  to  be  matched  to  achieve  optimum  performance  [1,2]. 

Consequently,  for  the  interpretation  of  radar  data  in  terms  of  a  changing  frequency 
content,  we  need  a  representation  of  our  data  as  a  function  of  both  time  and  frequency. 
The  purpose  of  this  report  is  therefore  to  stress  the  importance  of  alternative  methods 
which  have  been  given  little  attention  in  the  past,  namely,  the  joint  time-frequency 
representation  of  signals.  The  problem  of  detecting  accelerating  targets  was  looked  at 
from  two  computational  viewpoints:  1)  Correlator  method,  that  is,  de-chirping  followed 
by  Doppler  processing,  and  2)  Time-Frequency  methods.  The  correlator  method, 
although  computationally  costly,  is  optimum  (maximum-likelihood),  and  hence  serves 
as  a  benchmark.  Time-Frequency  methods  have  been  considered  before  in  a  study  done 
by  [4].  The  time-frequency  method  proposed  in  [4]  was  based  on  the  definition  of  a 
discrete-time  Wigner  distribution  given  in  [5].  This  method,  while  providing  insights 
into  the  problem,  has  the  following  drawbacks:  1)  it  is  sub-optimal,  2)  its  maximum 
unambiguous  range  of  measurable  velocities  is  half  of  that  of  the  correlator  method,  or 
of  the  conventional  Doppler  processing  method  \  3)  it  relies  on  an  exhaustive  search  of 
the  time-frequency  plane  in  much  the  same  way  as  the  correlator  method  relies  on  an 
exhaustive  search  of  the  velocity-acceleration  plane,  and  hence  could  not  be 
implemented  in  real-time. 

In  this  report,  we  derived  a  discrete-time  Wigner  time- frequency  distribution  which 
overcomes  drawbacks  1)  and  2).  This  distribution  was  independently  derived  by  one  of 

applicable  when  the  target  is  non-accelerating  or  mildly  accelerating,  see  [1]  for  more  details 
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the  authors  [1].  But  it  had  already  been  derived  by  Chan  [6]  in  an  effort  to  solve  the 
problem  of  aliasing.  Nevertheless,  the  optimality  properties  of  this  discrete-time 
Wigner  distribution  for  signal  detection  were  not  observed  in  [6].  Therefore,  in  the 
context  of  signal  detection,  this  discrete-time  Wigner  distribution  seems  new. 

The  performance  of  this  new  method  depends  on  the  signal-to-noise  ratio  and  the 
length  of  the  observation  interval  in  exactly  the  same  way  as  the  performance  of  the 
conventional  Doppler  processing  method  for  non-accelerating  targets  depends  on  the 
signal-to-noise  ratio  and  the  length  of  the  observation  interval.  Performance  analyses  of 
the  conventional  Doppler  processing  method  for  non-accelerating  targets  is  available  in 
the  literature,  for  example  [7].  Therefore,  the  performance  of  our  new  Wigner 
distribution  based  method  has  been  automatically  quantified.  This  new  Wigner 
distribution  has  led  to  the  conception  of  an  algorithm,  which  although  is  not  optimal, 
could  be  implemented  in  real-time  and  hence  would  overcome  drawback  3).  The  work 
presented  in  this  report  was  carried  out  as  part  of  a  collaborative  research  contract  from 
DREO  to  Raytheon  Canada  [1]. 

In  the  following  section  we  define  the  signal  model  on  which  the  algorithms  are  based. 
Section  3  discusses  the  maximum-likelihood  method  and  its  implementation  as  a 
correlator  or  as  dechirping  followed  by  Doppler  processing.  Section  4  discusses  three 
Wigner  distribution  based  methods,  one  of  which  is  shown  to  be  equivalent  to  the 
maximum-likelihood  method.  Section  5  identifies  an  algorithm  that  has  real-time 
feasibility  for  further  study.  Finally,  general  conclusions  are  summarized  in  Section  6. 
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2.  Signal  Model  and  Detection  and  Parameter 

Estimation  Problems _ 

In  this  study,  the  radar  is  assumed  to  be  of  the  Pulse  Doppler  type,  that  is,  it  sends  out  a 
uniform  train  of  RF  pulses  and  phase-coherently  receives  their  returns.  It  is  also 
assumed  that  the  receiver  has  pulse  compression  capability  so  that  a  pulse  return  from 
an  isolated  target  can  be  represented  by  a  single  sample  of  the  compressed  pulse  2. 

2.1  Single-Target  Scenario 

Suppose  the  radar  sends  out  N  pulses,  one  every  T  seconds,  and  there  is  a  target 
moving  with  constant  radial  acceleration.  Assume  that  the  range  walk  is  negligible,  that 
is,  the  change  in  range  during  the  obeservation  period  of  NT  seconds  is  negligible 
compared  to  the  radar  range  resolution  as  determined  by  the  width  of  the  compressed 
pulse.  Then  the  samples  of  the  N  range-compressed  pulses  taken  at  the  range  of  the 
targe  have  the  form 


(1)  r(n)  =  s(n)  +  v(n),  for  n  =  0, 1,2,3, . . . ,  (N  -  1), 

where 

(2)  s(n)  =  ae>{bo+bin+2b2n2\ 

is  the  noise-free  signal,  and  v(n)  is  a  sequence  of  independently  and  identically 
distributed  samples  of  complex  Gaussian  noise  with  zero  mean  and  variance  <r2. 

The  signal  parameters  a  and  &o  are  the  target  amplitude  and  phase  respectively.  The 
signal  parameters  b\  and  are  the  normalized  initial  radial  velocity  and  the  normalized 
radial  acceleration  respectively. 

2.1 .1  Normalized  Initial  Velocity  and  Acceleration 

The  normalized  initial  radial  velocity  is  defined  as 

...  .  ( 4nT\ 

where  u  is  the  initial  radial  velocity  in  meters/sec,  T  is  pulse  repetition 
interval  in  secs,  and  A  is  the  carrier  wavelength  in  meters.  Similarly,  the 
normalized  radial  acceleration  is  defined  as 

Alternatively,  the  transmitted  pulse  may  be  narrow  enough  to  achieve  the  same  effect. 
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(4) 


where  /  is  the  radial  acceleration  in  meters/sec/sec  towards  the  radar.  It  can  be 
seen  that  both  of  the  above  normalized  quantities  are  non-dimensional. 

To  obtain  these  relations  note  that  the  increment  in  range  d ,  measured  from 
the  beginning  of  the  observation  interval,  as  a  function  of  continuous-time  t ,  is 
d(t )  =  ut  +  l  ft2,  and  therefore  the  phase  increment  of  a  pulse  as  a  function 
of  discrete-time,  or  pulse  index,  is  given  by 


In  the  rest  of  the  discussion,  the  term  velocity  and  acceleration  refer  to 
normalized  radial  velocity  and  normalized  radial  acceleration 
respectively. 

2.1 .2  Per-Pulse  Signal-to-Noise  Ratio 

The  Per-Pulse  Signal-to-Noise  Ratio  is  defined  as 


(7)  SNRpulse=  02. 

The  performance  measures  of  any  detection  method  will  depend  on  SNRpuige 
and  the  number  of  pulses  N.  The  values  of  the  signal  parameters  bi  and  &2 
may  also  influence  performance.  We  are  interested  in  detection  methods 
whose  performance  will  increase  as  N  increases  regardless  of  the  values  of  bi 
and  &2- 


2.2  Multiple-Target  Scenario  -  Targets  Separated  in  Range 

When  there  are  multiple  targets,  as  long  as  they  are  separated  from  one  another  in  range 
by  at  least  the  radar  range  resolution  as  determined  by  the  width  of  the  compressed 
pulse,  a  model  of  the  above  form  will  be  valid  for  each  of  them. 

Thus,  if  there  are  two  targets  separated  in  range,  the  samples  of  the  N 
range-compressed  pulses  taken  at  the  corresponding  ranges  of  the  targets  have  the  form 
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(8)  n(n)  =  $i(n)  +  vi  (n),  and 

(9)  r2(n)  =  s2(n)  +  v2(n) 


for  n  —  0, 1,2,3, . . . ,  (iV  —  1), 


where  sx(n)  and  s2 (n)  are  the  noise-free  signals,  and  i>i(n)  and  v2(n)  are  each 
sequences  of  independently  and  identically  distributed  (iid)  samples  of  complex 
Gaussian  noise  with  zero  mean  and  variance  a2. 

Because  of  the  independence  of  Ui(n)  and  v2(n),  the  model  decouples  into  two  single 
target  models. 


2.3  Multiple-Target  Scenario  -  Targets  Coincident  in  Range 

When  there  are  multiple  targets  coincident  in  range,  the  samples  of  the  N 
range-compressed  pulses  taken  at  the  range  of  the  targets  have  the  form 


(10)  r(n )  =  s(n)  +v(n ), 


for  n  =  0, 1,2,3,. . (iV  —  1), 


where  s(n)  is  the  combined  noise-free  signal  due  to  all  targets  and  v(n)  is  a  sequence 
of  independently  and  identically  distributed  (iid)  samples  of  complex  Gaussian  noise 
with  mean  zero  and  variance  a2.  If  there  are  K  targets,  the  combined  noise-free  signal 
has  the  form 


(ID 


K 

s(n )  =  y^afce^M+6Mn+K2n2), 

k—  1 


where  and  o  are  the  amplitude  and  phase  respectively  of  the  kfi1  target,  and  b^x 
and  2  are  the  normalized  initial  radial  velocity  and  the  normalized  radial  acceleration 
respectively  of  the  kfi  target. 


2.4  Target  Amplitude  and  Phase  Models 

The  target  phase  bo  is  assumed  to  be  a  random  variable  uniformly  distributed  over 
[— 7r,  7r].  The  target  amplitude  a  may  be  assumed  to  be  a  constant  (Marcum,  or 
Swerling-0)  or  a  random  variable  [7].  Commonly  assumed  probability  functions  are  the 
Rayleigh  density  (Swerling-1)  and  Rice  density. 
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2.5  Detection  Problem 


Given  the  samples  of  the  N  range-compressed  pulses  taken  at  a  particular  range,  the 
detection  problem  is  to  choose  between  the  following  hypotheses: 

•  Noise-only  Hypothesis  -  Hq,  that  is,  the  received  signal  r(n)  is  due  to  noise  alone, 

(12)  r(n)  —  v(n ), 

•  Signal-and-Noise  Hypothesis  -  Hi,  that  is,  the  received  signal  r(n)  is  due  to  signal 
and  noise, 

(13)  r(n)  —  s(n)  +  t?(n), 

where  the  signal  s(n)  may  be  due  to  the  presence  of  one  or  more  targets  at  the  range  of 
interest. 

When  there  is  more  than  one  target,  the  ‘Signal-and  Noise’  hypothesis  can  be 
subdivided  on  the  basis  of  the  number  of  hypothesized  targets,  or  chirp  signals,  as  ‘  1 
target  and  noise’,  ‘2  targets  and  noise’,  ‘3  targets  and  noise’,  and  so  on. 

2.6  Parameter  Estimation  Problem 

For  the  purpose  of  accurately  tracking  the  accelerating  targets  we  need  to  estimate  their 
parameters,  especially  their  kinematic  parameters  such  as  initial  velocity,  median 
velocity,  and  acceleration.  Therefore,  the  problems  of  detecting  the  targets  and 
estimating  their  parameters  must  be  posed  as  a  combined  problem. 
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3.  The  Maximum  Likelihood,  or  Generalized  Likelihood 
Ratio  Test,  Method 


In  this  section,  the  maximum  likelihood,  or  generalized  likelihood  ratio  test,  method  of 
detecting  an  accelerating  target  and  estimating  its  parameters  is  described.  This  method 
is  asymptotically  optimum,  and  hence  will  sometimes  be  referred  to  as  the  optimum 
method  in  this  report  [8-12]. 

Recall  from  Section  2  that  the  detection  problem  is  to  choose  between  the  Noise-Only 
Hypothesis  -  H0  and  the  Signal-and-Noise  Hypothesis  -  Hi. 

Define 


N-l  N- 1 

(14)  A (d,  co,  ci,  C2)  =  lr(n)|2  ~~  ^2  |r(n)  —  de?(C0+Cin+2C2n2) 

n=0  n=0 

Then 


(15) 


Hi 

max  ...  .  > 

j  „  .  „  A(d,flo,Ci,p2) 

a,co,ci,C2  < 

H0 


for  some  threshold  7?,  is  the  maximum  likelihood  (ML),  or  generalized  likelihood  ratio 
test  (GLRT),  method.  The  maximizing  values  of  d ,  Co,  Ci,  and  C2  are  the  maximum 
likelihood  estimates  of  a,  bo,  and  62  (see  Section  2). 

It  can  be  shown  that  the  above  test  can  take  the  equivalent  form 


(16) 


Hi 


max 

Cl,C2 


A(ci,c2)  ^ 

Ho 


7, 


where 


(17) 


A(ci,c2) 


N- 1 

^2  r{n)e~i(c'n+ iC2"2) 

n=0 
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and  7  is  a  threshold  (which  can  be  related  to  rj).  Here  the  maximizing  values  of  ci  and 
C2  are  the  maximum  likelihood  estimates  of  &i  and  62  respectively  [1],  The  optimum  d 
and  Co  are  given  in  terms  of  the  optimum  c\  and  C2  by  [1  ] 

1  W_1 

(1 8)  <feJC°  =  r(n)e~j{c'nh*rW)- 

n= 0 

Thus  the  ML,  or  GLRT,  method  can  be  implemented  as  a  correlator  in  the  initial 

velocity  -  acceleration  plane.  This  method  is  also  known  as  the  dechirp-Doppler 

■  1  2 

method,  i.e.,  first  multiplying  by  e~J2C2n  to  obtain  a  pure  complex  exponential,  or 
nearly  so,  so  that  there  will  be  little  or  no  loss  due  to  Doppler  spreading  when 
estimating  the  frequency  of  the  pure  complex  exponential  by  the  conventional  Doppler 
processing  method. 

The  notion  of  output  SNR  of  the  method  comes  from  the  observation  that  under  the 
‘noise-only’  hypothesis  A(ci,C2)  is  a  central  Chi-Squared  random  variable  of  two 
degrees  of  freedom,  for  all  C\  and  C2,  whereas  under  the  ‘signal-and-noise’  hypothesis, 
-^2  A(&i,  62)  is  a  non-central  Chi-Squared  random  variable  of  two  degrees  of  freedom 

and  non-centrality  parameter  2 N  (^ )  [6,7]. 

We,  therefore,  consider 

(19)  SNRout  —  N  • 

as  this  can  be  used  to  approximately  compute  the  probability  of  detection,  which  is 
sometimes  assuming  that  the  maximum  of  A(ei,  c?)  occurs  at  c±  —  b\  and  c<i  —  f>2- 
Note  that  SNRoUt  is  a  monotonically  increasing  function  of  N,  but  more  specifically,  it 
is  a  linearly  increasing  function  of  N. 

3.1  The  Cost  of  Computation  for  an  Exhaustive  Search 

The  cost  of  computing  A(ci,  C2),  by  brute-force,  over  K  values  of  C2  and  M  values  of 
ci  is 

(20)  KMO  (N) . 

It  is  clear  that  for  a  given  C2,  we  can  compute  A(ci,  C2)  for  M  uniformly  spaced 
discrete  values  of  C\  in  the  range  [— 7r,  7t]  by  using  the  FFT  algorithm,  provided 
M  >N.  Thus,  if  we  desire  to  calculate  A(ci,  C2)  over  K  values  of  C2  and  M  values  of 
Ci,  uniformly  spaced  over  [— 7r,  tt]  then  the  computational  cost  is 
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(21) 


KO  (M  log  M) . 
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4.  A  Wigner-Distribution  Formulation  of  the  Optimum 
Detection  Problem  for  a  Discrete-Time  Chirp  Signal 


In  this  section  we  formulate  the  maximum-likelihood,  or  optimum,  detection  of  a 
discrete-time  chirp  signal  in  additive  white  Gaussian  noise  in  terms  of  a  Wigner 
distribution  of  the  received  signal  r(n).  We  do  this  by  showing  that  Ar(ci,  c2)  can  be 
interpreted  as  the  integral  along  a  straight  line  on  the  time-frequency  plane  of  an 
appropriately  defined  Wigner  distribution  of  the  discrete-time  signal  r(n). 

Thus  the  maximization  over  c\  and  c2  is  transformed  into  maximization  over  straight 
lines  whose  intercept  and  slope  parameters  correspond  to  ci  and  c2  respectively.  At  first 
this  does  not  seem  to  offer  any  advantages  over  direct  maximization  of  Ar(ci,  c2). 
However,  the  Wigner  distribution  may  suggest  non-exhaustive,  fast  methods  of 
searching  and  identifying  the  region  in  the  time-frequency  plane  where  the  target’s 
signal  exists.  This  is  especially  so  because,  in  the  absence  of  noise  or  any  other 
interference,  the  Wigner  distribution  of  a  discrete-time  chirp  signal  is  concentrated  on  a 
straight  line  of  the  time-frequency  plane.  A  method  based  on  this  observation  is 
identified  in  Section  5  for  further  investigation. 

In  the  presence  of  interference  of  kinds  other  than  white  Gaussian  noise,  the  Wigner 
distribution  may  suggest  sub-optimal,  robust  methods  of  detection.  For  example,  if 
some  unwanted  interference  is  known  to  be  concentrated  in  a  certain  region  of  the 
time-frequency  plane,  then  that  interference  can  be  filtered  out  and  a  cleaner  signal  can 
be  synthesized  for  further  processing  by  more  conventional  means. 

4.1  Three  Definitions  of  Wigner  Distribution  for  Discrete-Time 
Signals 

Given  a  discrete-time  signal  r(n),  we  consider  three  definitions  of  auto  time-frequency 
distributions  which  have  some  properties  analogous  to  those  of  the  classical  Wigner 
distribution  defined  for  continuous-time  signals  [13]. 

It  is  straightforward  to  verify  that  these  are  real-valued  functions. 

4.1.1  Type-1  Wigner  Distribution 

The  type-I  Wigner  distribution  W* (n,  8)  is  defined  as 


(22)  Wf(n?8)  =  ^^r(n-\-k)r*(n  —  k)e 

k 

where  n  is  integer- valued  and  6  is  real-valued.  Note  that  this  is  the  same  as  the 
definition  provided  in  [5]  except  for  the  missing  scaling  factor  2  at  the  front. 
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It  is  easy  to  verify  that  the  type-I  Wigner  distribution  is  a  periodic  function  of 
9  with  period  7r.  This  property  will  be  discussed  more  in  the  following 
subsections. 

For  a  signal  r(n)  that  is  zero  outside  0  <  n  <  (N  —  1),  the  type-I  Wigner 
distribution  is  zero  outside  0  <  n  <  (N  -  1). 

4.1 .2  Type-11  Wigner  Distribution 

The  type-II  Wigner  distribution  W//(n,  9)  is  defined  as 


(23)  W^(n,  0)  =  ^  r(n  +  k  +  l)r*(n  -  k)e~^2k+1^, 

k 

where  n  is  integer- valued  and  9  is  real- valued  [1]. 

It  is  easy  to  verify  that  the  type-II  Wigner  distribution  is  a  periodic  function  of 
9  with  period  2n.  It  is  also  easy  to  verify  that 

(24)  W^in,  9  +  tt)  -  -W^(n,  9 ). 

These  properties  will  be  discussed  more  in  the  following  subsections. 

For  a  signal  r(n)  that  is  zero  outside  0  <  n  <  (N  —  1),  the  type-II  Wigner 
distribution  is  zero  outside  0  <  n  <  (N  —  2). 

4.1 .3  Type-lll  Wigner  Distribution 

The  type-Itl  Wigner  distribution  Wfn(n ,  9)  is  defined  in  terms  of  the  type-I 
and  type-II  Wigner  distributions  as  follows  [1]. 


(25) 


W^n/2,0) 

W”({n-l)/2t9) 


for  even  n> 
for  odd  n. 


It  is  easy  to  verify  that  the  type-III  Wigner  distribution  is  a  periodic  function 
of  9  with  period  2n.  This  property  will  be  discussed  more  in  in  the  following 
subsections. 

For  a  signal  r(m)  that  is  zero  outside  0  <  m  <  (N  —  1),  the  the  type-III 
Wigner  distribution  is  zero  outside  0  <  n  <  2(JV  —  1).  However,  if  we 
consider  even  n  to  correspond  to  integer  values  of  time  and  odd  n  to 
correspond  to  half-integer  values  of  time  then  the  type-HI  Wigner  distribution 
is  zero  outside  the  time  range  0  <  m  <  (N  —  1). 
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4.2  Computing  Ar(ci,  C2)  via  the  Wigner  Distributions 

In  Appendix  A  we  have  shown  that 


N- 1  N~  2 

(26)  Ar (*,<*)  =  £  Wrf («, ci  +  c2n)  +  ^  WrJ/ (n, ct  +  - c2  +  c2n) , 

n=0  n=0 

N- 1  AT— 2 

(27)  Ar(ci,  c2)  =  ^  Wjn(2n,  cx  +  - c22n )  +  ^  Wr777(2n  + 1,  ci  +  -c2(2n+ 1)), 

n=0  n=0 

2N—2  . 

(28)  Ar(ci,c2)  =  ^2  M/r///(n-Ci  +  2C2n)- 

77—0  ** 

Thus  we  can  calculate  Ar(ci, C2)  by  taking  the  type-111  Wigner  distribution  Wj.n{n, 6) 
of  the  discrete-time  received  signal  r(n)  and  integrating  it  along  the  line  with  intercept 
c\  (value  of  6  at  n  =  0)  and  slope  ~c2  (increment  in  6  per  unit  increment  in  n). 

This  implies  that  we  can  perform  optimum  detection  by  the  rule 


(29) 


max 

ci,c2 


2N—2 


Ew‘ 


n— 0 


111  (n,  a  +  -c2n) 


) 


Hi 

> 

< 

Ho 


7? 


for  some  appropriately  chosen  threshold  7,  where  Ho  and  Hi  are  the  ‘noise  only’  and 
‘signal  plus  noise’  hypotheses  defined  in  Section  2. 

The  type-III  Wigner  distribution  was  independently  derived  by  one  of  the  authors  [1]. 
But  it  had  already  been  derived  by  Chan  [6]  in  an  effort  to  solve  the  problem  of 
aliasing.  Nevertheless,  the  optimality  properties  of  this  discrete-time  Wigner 
distribution  for  signal  detection  were  not  observed  in  [6].  Therefore,  in  the  context  of 
signal  detection,  this  discrete-time  Wigner  distribution  seems  new. 


4.3  Three  Definitions  of  Radon-Wigner  Transform  and  Detection 
Methods  Based  on  them  for  Discrete-Time  Chirp  Signals 

Based  on  the  results  of  Section  4.2,  we  can  define  the  Radon-Wigner  transform  of  a 
discrete-time  signal  in  three  ways  and  conceive  of  three  detection  methods  for 
discrete-time  chirp  signals.  In  the  literature  the  term  ‘Hough’  transform  is  sometimes 
used  to  refer  to  what  we  have  called  ‘Radon’  transform.  A  more  detailed  description  of 
the  Radon  transform  can  be  found  in  f  14]. 
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4.3.1  Type-1  Radon-Wigner  transform  and  a  detection  method  based  on  it 


We  shall  refer  to 


(30) 


N- 1 

RW* (*,<*)  =  ^2  W/(n,ci  +c2n) 

n=0 


as  the  Type-I  Radon-Wigner  transform  of  r(n )  [1], 

We  may  use  the  rule 

Hi 

<31)  ( Zl  <  '''• 

Ho 

for  some  appropriately  chosen  threshold  77,  to  detect  discrete-time  chirp 
signals. 

Equation  (30)  can  be  summarized  as  follows.  To  optimally  detect  an 
accelerating  target  in  white  Guassian  noise  one  should  compute  the  wigner 
distribution  of  the  observed  data  along  all  lines  in  the  time-frequency  plane. 
The  maximum  integrated  value  is  then  compared  to  a  threshold  and  a 
detection  is  declared  if  the  threshold  is  exceeded.  It  is  interesting  to  note  that 
this  transformation  of  a  signal  via  line  integration  to  another  2-D  signal  is  a 
Radon  transform  [14].  In  practice,  one  usually  has  some  a  priori  knowledge 
about  the  range  of  possible  instantaneous  frequency  values  which  allows  one 
to  restrict  the  line  integration  to  specific  regions  of  the  time-frequency  plane. 

4.3.2  Type-ll  Radon-Wigner  transform  and  a  detection  method  based  on 
it 

We  shall  refer  to 


(32)  RW”{cu c2)  =  J2  ci  +  \<*  +  c*n) 

n— 0 

as  the  Type-II  Radon-Wigner  transform  of  r(n)  [1]. 

We  may  use  the  rule 


DREO  TM  2001-141 


13 


(33) 


(max 
ci,c2 


Hx 

> 

< 

Ho 


for  some  appropriately  chosen  threshold  7/J,  to  detect  discrete-time  chirp 
signals. 

4.3.3  Type-Ill  Radon-Wigner  transform  and  a  detection  method  based  on 
it 

We  shall  refer  to 


2N—2 

(34)  RW?II(cl,c2)  =  Wln(n,ci  +  -c2n) 

n=0 

as  the  Type-III  Radon-Wigner  transform  of  r(n)  [1]. 

A  rule  based  on  this  to  detect  discrete-time  chirp  signals  has  been  stated  in 
Section  3.  We  restate  it  as 


°5)  ( ““ 

Ho 

4.3.4  A  Common  Error  and  its  Consequences 

It  is  a  common  error  to  assume 


(36) 

(37) 


JV-l 

Ar(ci,c2)  =  2  23  W/(n>ci  +  c2 n), 


n=0 

=  2  RWl(cuc2), 


and  claim  that  optimum  detection  can  be  performed  by  the  rule 
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(38) 


(  “  ***}&*>)  <  T 

Ho 

This  error  was  made  in  [4],  An  historical  reason  for  how  this  may  have 
happened  is  given  in  Appendix  B. 

This  error  has  the  following  negative  consequences: 

•  A  50%  reduction  in  the  range  of  velocities  that  can  be  unambiguously 
measured, 

•  A  loss  in  the  asymptotic  performance  as  measured  by  the  output 
signal-to-ratio  relative  to  the  optimum  method. 


4.3.5  The  Range  of  Unambiguously  Measurable  Velocities  of  the  Type-1 
and  Type-Ill  Wigner  Distribution  Based  Methods 

As  was  shown  in  Section  4.1.1,  the  period  of  the  type-I  Wigner  distribution  is 
7r.  Therefore  the  signals 


(39) 

(40) 


aej(fco+6in+l62n2);  wd 
a€j(bo-\-(bi-\-7r)n-\-lb2Ti2)  ^ 


will  give  rise  to  the  same  type-I  Wigner  distribution.  In  otherwords,  if  we 
detect,  using  the  type-I  Wigner  distribution,  the  presence  of  a  signal  in  the 
time-frequency  plane  with  intercept  and  slope  parameters  c\  and  C2 
respectively,  then  we  cannot  resolve  the  ambiguity  as  to  whether  the  physical 
signal  was 


aeJ'(co+cin+2C2f,2>,  or 
ae?(C0+(ci+*)*»+2  C2fl2) 

without  external  information. 

In  contrast,  as  was  shown  in  Section  4.1.3,  the  period  of  the  type-III  Wigner 
distribution  is  2n.  Therefore  the  signals 


(41) 

(42) 
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(43) 

(44) 


Oei(«o+6in+ifc,»a)j  ^ 

ae?  (&0 +  (*>1 +*)™+  \  bin2  ) 


will  give  rise  to  different  type-III  Wigner  distributions. 

If  we  detect,  using  the  type-III  Wigner  distribution,  the  presence  of  a  signal  in 
the  time-frequency  plane  with  intercept  and  slope  parameters  c\  and  c<i 
respectively,  then  there  is  no  ambiguity  as  to  whether  the  physical  signal  was 


(45)  ^(co+cm+icjn2),  or 

(46)  aei(co + (ci +7r)n+  \  c2n2 )  ^ 

But  we  would  of  course  need  external  information  to  resolve  the  ambiguity  as 
to  whether  the  physical  signal  was 


aeJ(Co+Cl"+2C2"2),  or 
aei(co + (ci  +2ir  )n+  \  c2  n2 ) 

4.3.6  The  Output  Signal-to-Noise  Ratio  of  the  Type-1  and  Type-III  Wigner 
Distribution  Based  Methods 

Owing  to  the  equivalence  shown  in  Section  4.2,  the  output  Signal-to-Noise 
Ratio  (SNR)  of  the  type-III  Wigner  distribution  based  method  is  the  same  as 
that  of  the  optimum  maximum  likelihood  method  described  in  Section  3.  The 
peformance  results  on  the  latter  is  available  in  the  literature. 

It  was  claimed  in  [15]  that  the  output  Signal-to-Noise  Ratio  (SNR)  of  the 
type-I  Wigner  distribution  based  method  was  less  by  3dB  than  that  of  the 
N(^)2  (i.e.,  type-III  Wigner  distribution  based  method). 

4.4  The  Wigner  Distributions  of  the  Discrete-Time  Chirp  Signal 

In  view  of  the  results  of  Section  4.2,  it  is  of  interest  to  derive  the  type-I,  type-II,  and 
type-III  Wigner  distributions  for  the  discrete-time  chirp  signal 


(47) 

(48) 


(49)  s(n)  =  aeJ'(&0+6l”+262n2),  for  n  =  0, 1, . . . ,  (N  -  1), 

defined  in  Section  2.  These  derivations  have  been  done  in  Appendix  C  assuming  a  =  1. 
The  results  stated  here  are  obtained  by  scaling  those  derived  in  the  appendix  by  a2. 
Note  that  the  Wigner  distributions  do  not  depend  on  bo- 
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4.4.1 


The  Type-1  Wigner  Distribution  of  s(n ) 

The  Type-I  Wigner  Distribution  Wl (n,  9)  of  s(n)  is  given  by  [1] 


Wl{n,9) 

(50) 


„  f  2  min(n,  N  —  1  —  n)  + 1  if  0  =  &i  +  &2«  mod  it, 

a  |  otherwjse 


for  0  <n<(N  —  1),  and  it  is  zero  outside  this  range. 


In  a  3-dimensional  plot,  it  has  a  ridge  along  the  lines  given  by 


(51)  0  =  bi  +  b^n  mod  n. 

The  height  of  the  ridge  is  given  by 


(52)  2min(n,iV  —  1  —  n)  +  1. 

4.4.2  The  Type-ll  Wigner  distribution  of  s(n ) 

The  Type-II  Wigner  Distribution  W/1  (n,  9)  of  s(n)  is  given  by  [1] 


W3n(n,9)  =  a2 
(53) 


2  min(n,  N  —  2  —  n)  +  2 
—(2  min  (n,N  —  2  —  n)  +  2) 

8in[(g—  (hi  +  §&2+i>2rc))(2  min(n,JV— 2— n)+2)] 
sin[^-(6i+^62+&2n)] 


if  9  =  fa  +  ^62  +  i>2 n  mod  27r, 
if  9  =  7r  +  61  +  mod  2n, 

otherwise. 


for  0  <  n  <  (N  —  2),  and  it  is  zero  outside  this  range. 


In  a  3-dimensional  plot,  it  has  a  ridge  along  the  lines  given  by  [1] 


(54)  9  =  bi  -f  mod  27r. 

The  height  of  the  ridge  is  given  by 


(55)  2min(n,iV  —  2  —  n)  +  2. 

In  a  3-dimensional  plot,  it  also  has  a  valley  along  the  lines  given  by 
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(56) 


0  —  7r  +  fa  +  -fa  +  fan  mod  n. 


The  depth  of  the  valley  is  given  by 


(57)  2  min(n,  JV  —  2  —  n)  +  2. 

4.4.3  The  Type-Ill  Wigner  distribution  of  s(n ) 

The  Type-Ill  Wigner  Distribution  W//J(n,  0)  of  s(n )  is  given  by 


WgH(n,  6)  =  a 2 
(58) 


2min(|,JV  — l-f)  +  1 
(-1)”  [2min(s,JV-l-  f )  +  l] 

sin[(fl— (6i  +  |t>2Ti))(2min(^,jV— 1— If  )+l)j 

sin[fl-(6i  +  |62n)] 


if  0  =  bi  +  \b2r1  mod  27r, 
if  6  =  7r  +  bi  +  fan  mod  27r, 

otherwise. 


for  0  <  n  <  2 (N  —  1),  and  it  is  zero  outside  this  range. 


In  a  3-dimensional  plot,  it  has  a  ridge  along  the  lines  given  by 


(59) 


0  =  6i  +  ~fan  m0(^  27t. 


The  height  of  the  ridge  is  given  by 


(60)  min(n,2(Ar  —  1)  —  n)  +  1. 

In  a  3-dimensional  plot,  it  also  has  an  oscillatory  structure  along  the  lines  give 
by 


(61)  6  =  7r  -I-  61  +  -zb^n  mod  2n. 

The  period  of  the  oscillation  is  one  unit  and  the  amplitude  of  the  oscillation  is 
given  by 


(62) 


min(n,2(iV  —  1)  —  n)  +  1. 
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4.5  A  Numerical  Example 


The  example  of  a  two- target  scenario,  observed  with  N  =  128,  described  in  Section  2, 
is  studied  here  using  the  three  types  of  Wigner  distribution  based  methods.  The 
normalized  initial  velocity  and  normalized  acceleration  for  target  1  are  0  and  0, 
respectively.  For  the  target  2,  the  normalized  initial  velocity  and  normalized 
acceleration  are  1  and  0.01,  respectively. 

Figure  1  shows  the  contour  plot  for  the  type-I  Wigner  distribution  of  the  noise-free 
signal  whereas  Figure  2  shows  the  surface  plot  of  the  noise-free  signal.  Note  that  the 
time-frequency  of  the  accelerating  target’s  Wigner  distribution  is  concentrated  along  a 
linear  segment  in  the  time-frequency  plane.  This  linear  segment  corresponds  to  the 
instantaneous  frequency.  Figure  3  shows  the  contour  plot  for  the  type-I  Wigner 
distribution  of  the  noisy  signal  and  Figure  4  shows  the  corresponding  surface  plot.  For 
the  Wigner  distribution  plotted  in  Figures  3  and  4,  the  targets  were  contaminated  with 
the  random  white  Guasian  noise,  such  that  the  SNR  was  5  dB.  It  is  apparent  that  the 
type-I  Wigner  distribution  is  a  periodic  function  with  period  7r.  Figures  5  and  6  show 
the  type-I  Radon- Wigner  transform  of  the  noise-free  signal  showing  target  1  and  target 
2,  respectively.  Figures  7  and  8  show  the  type-I  Radon- Wigner  transform  of  the  noisy 
signal  showing  target  1  and  target  2,  respectively.  The  plots  of  the  Wigner  distriburion 
and  Radon- Wigner  transform  show  that  an  accelereting  target  can  be  detected  as  easily 
as  a  non-accelerating  target. 

Figure  9  shows  the  contour  plot  for  the  type-II  Wigner  distribution  of  the  noise-free 
signal  and  Figure  10  shows  the  surface  plot  of  the  noise-free  signal.  Figure  1 1  shows 
the  contour  plot  for  the  type-II  Wigner  distribution  of  the  noisy  signal  whereas  Figure 
12  shows  the  corresponding  surface  plot.  It  is  evident  that  the  type-II  Wigner 
distribution  is  a  periodic  function  with  period  2n  and  Wf1  (n,  0  +  7r)  =  — W/7(n,  0) 
(see  Section  4.1.2).  Figures  13  and  14  show  the  type-II  Radon- Wigner  transform  of  the 
noise-free  signal  showing  target  1  and  target  2,  respectively.  Figures  15  and  16  show 
the  type-II  Radon- Wigner  transform  of  the  noisy  signal  showing  target  1  and  target  2, 
respectively.  The  plots  of  the  Wigner  distribution  and  the  Radon- Wigner  transforms 
show  that  an  accelerating  target  can  be  detected  as  easily  as  a  non- accelerating  target. 

Figure  17  shows  the  contour  plot  for  the  type-III  Wigner  distribution  of  the  noise-free 
signal  and  Figure  18  shows  the  contour  plot  for  the  type-III  Wigner  distribution  of  the 
noisy  signal.  It  is  apparent  that  the  type-III  Wigner  distribution  is  a  periodic  function 
with  period  27r  and  this  property  is  discussed  more  in  Section  4.3.4.  Figures  19  and  20 
show  the  type-III  Radon- Wigner  transform  of  the  noise-free  signal  showing  target  1  and 
target  2,  respectively.  Figures  21  and  22  show  the  type-in  Radon- Wigner  transform  of 
the  noisy  signal  showing  target  1  and  target  2,  respectively.  The  plots  of  the  Wigner 
distribution  and  the  Radon- Wigner  transforms  show  that  an  accelerating  target  can  be 
detected  as  easily  as  a  non-accelerating  target. 
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4.6  Type-1  Wigner  Distribution  and  its  Radon  Transform 

4.6.1  Wigner  Distribution 


NoiseFreeWDI 


Discrete  Time 


Figure  1:  The  Type-1  Wigner  Distribution  of  the  Noise-Free  Signal. 


NoiseFreeWDI 


Figure  2:  The  Type-1  Wigner  Distribution  of  the  Noise- Free  Signal. 


20 


DREO  TM  2001-141 


NoisyWDI,  SNR  =  5  dB 


0  20  40  60  80  100  120 

Discrete  Time 


Figure  3:  The  Type- 1  Wigner  Distribution  of  the  Noisy  Signal. 


NoisyWDI,  SNR  =  5  dB 


Figure  4:  The  Type- 1  Wigner  Distribution  of  the  Noisy  Signal. 


4.6.2  Radon  Transform 

NoiseFreeRadonl,  Target  1 
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Figure  5:  The  Type- 1  Radon-Wigner  Transform  of  the  Noise-Free  Signal  Showing  Target  1. 


NoiseFreeRadonl,  Target  2 
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Normalized  Acceleration 


Figure  6:  The  Type-1  Radon-Wigner  Transform  of  the  Noise-Free  Signal  Showing  Target  2. 
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Figure  7:  The  Type-1  Radon-Wigner  Transform  of  the  Noisy  Signal  Showing  Target  1. 


NoisyRadonl,  Target  2,  SNR  =  5  dB 
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Figure  8:  The  Type-1  Radon-Wigner  Transform  of  the  Noisy  Signal  Showing  Target  2. 
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4.7  Type-ll  Wigner  Distribution  and  its  Radon  Transform 

4.7.1  Wigner  Distribution 


20  40  #0  60  100  120 

Discrete  Time 


Figure  9:  The  Type-ll  Wigner  Distribution  of  the  Noise-Free  Signal. 


Figure  10:  The  Type-ll  Wigner  Distribution  of  the  Noise-Free  Signal. 
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4.7.2 


Radon  Transform 


Figure  13:  The  Type-H  Radon-Wigner  Transform  of  the  Noise-Free  Signal  Showing  Target  1. 


Figure  14:  The  Type-1 1  Radon-Wigner  Transform  of  the  Noise-Free  Signal  Showing  Target  2. 
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NoisyRadon2,  Target  1 ,  SNR  =  5  dB 
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Figure  15:  The  Type-ll  Radon-Wigner  Transform  of  the  Noisy  Signal  Showing  Target  1. 


NoisyRadon2,  Target  2,  SNR  =  5  dB 


0.0075  o.ooe  o.ooes  0.009  0.0095  0.01  0.0105  0.011  0.011s  0.012  0.0125 


Normalized  Acceleration 

Figure  16:  The  Type-ll  Radon-Wigner  Transform  of  the  Noisy  Signal  Showing  Target  2. 
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4.8  Type-Ill  Wigner  Distribution  and  its  Radon  Transform 

4.8.1  Wigner  Distribution 
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Discrete  Time 


Figure  17:  The  Type-Ill  Wigner  Distribution  of  the  Noise-Free  Signal. 
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Figure  18:  The  Type-IH  Wigner  Distribution  of  the  Noisy  Signal. 
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4.8.2  Radon  Transform 


Figure  19:  The  Type-ill  Radon-Wigner  Transform  of  the  Noise-Free  Signal  Showing  Target  1. 
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Figure  20:  The  Type- III  Radon-Wigner  Transform  of  the  Noise- Free  Signal  Showing  Target  2. 
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NoisyRadon3,  Target  1,  SNR  =  5  dB 


Figure  21:  The  Type -Hi  Radon-Wigner  Transform  of  the  Noisy  Signal  Showing  Target  1. 


Figure  22:  The  Type-Ill  Radon-Wigner  Transform  of  the  Noisy  Signal  Showing  Target  2. 
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5.  A  Detection  Method  with  Real-Time  Feasibility 


In  this  section  a  new  method  of  detecting  discrete-time  chirp  signals  and  estimating 
their  parameters  is  described.  It  should  be  feasible  to  implement  the  method  in 
real-time.  Therefore,  the  method  deserves  further  study. 

This  method  is  motivated  by  the  type-II  Wigner  distribution  defined  in  Section  4.  More 
specifically,  it  makes  use  of  the  type-II  Wigner  distribution  evaluated  only  at  the  center 
of  the  observation  interval.  The  cost  of  computation  can  be  shown  to  be  of  the  order  of 
that  of  a  Fourier  transform.  We  first  consider  detecting  a  single  chirp  signal.  This  will 
set  the  stage  for  considering  the  detection  of  multiple  chirp  signals. 


5.1  Detecting  a  Single  Chirp  Signal 

Recall  from  Appendix  A  that  for  the  discrete-time  chirp  signal 


(63) 


l  0 


ifO<n<  (AT  — 1), 
otherwise, 


the  type-II  Wigner  distribution  is  given  for  0  <  n  <  (N  —  2)  by 


wan(n,e) 

(64) 


'  2min(n,AT  — 2  — n) +2 
i  —  (2min(n,  N  —  2  —  n)  +  2) 

8in[(fl—  (6i+|62+&2n))(2min(n,JV— 2— w)+2)| 
k  sin[0-(6i-t-i&2+62n)] 


if  8  =  bi  +  |&2  +  62 n  mod  2tt, 
if  &  =  7r  +  61  +  \t>2  +  b^n  mod  27r, 

otherwise. 


For  a  fixed  n,  the  maximum  of  W^r(n,  6)  with  respect  to  6  occurs  when 


(65)  0  =  &i  +  +  b2n  mod  27r. 

Particularly,  when  N  is  even,  for  n  —  the  maximum  of  W/7(n,  0)  with  respect  to 
8  occurs  when 


(66)  8  =  h  +  h  mod  2?r- 

Hereafter  we  shall  assume  that  N  is  even. 3 
3For  N  odd  we  must  consider  the  type-II  Wigner  distribution  at  either  n  =  or  n  =  . 
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Therefore  in  the  presence  of  noise,  provided  it  is  not  too  large,  the  chirp  signal  can  be 
detected  by  the  rule 


(67) 


Hi 

> 

Ho 


for  some  appropriate  7m,  where  Hi  and  Hq  are  the  ‘signal-plus-noise’  and 
‘noise-only’  hypotheses.  If  a  signal  is  detected,  then  assuming  that 


(68) 


<  7T, 


the  maximizing  9 ,  say  Cm ,  can  be  taken  as  an  estimate  of  the  median  frequency 
bm  =  h  -f  &2  (^p)  of  the  chirp  signal  over  the  observation  interval 
n  “  0, . . . ,  (iV  —  1). 

If  an  estimate  of  the  frequency-rate  &2  is  desired,  it  can  be  obtained  by 


max 

(69)  ci,C2  Ar(ci,c2), 

ci  +  c2  (^“)  =  Cm 

that  is  by  maximizing  Ar(ei,  c2)  with  respect  to  C\  and  c2  under  the  constraint  that 
Ci  +  c2  —  Cm ,  where  Cm  is  the  estimate  of  the  median  frequency  obtained  based 

on  W”(^,0)  and 


(70) 


Ar(ci,c2)  = 


JV-l 


2  r^e 


-j(cin+\c2n2) 


n=0 


We  can  restate  the  latter  maximization  problem  as 


(71) 


max 

c2 


Ar(cm-C2(^y^),C2). 


For  convenience  we  denote 
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I 


(72) 

(73) 


r{Cm,c2)  —  A T  {^m  —  c2  ^  2  ^  > c2^  » 

=  |^1r(n)e-^((ft»^(1^1)),H-i«w*) 


n=0 


A  better  trade-off  between  the  probabilities  of  ‘miss’  and  ‘false  alarm’  may  be  obtained 
by  using  a  lower  threshold  7m  in  the  first  test  and  turning  the  second  maximization 
problem  into  another  test  as4 


Hi 

(74)  Mwa)  J  7a, 

H0 

for  some  appropriate  7a.  This  idea  turns  out  to  be  useful  for  detecting  multiple  chirp 
signals. 


5.2  Detecting  Multiple  Chirp  Signals 

The  Wigner  distribution  of  a  superposition  of  signals  contains  cross  terms,  or 
interferences,  between  the  signals.  Suppose  Si(n)  and  S2(n)  are  two  signals,  and 
s(n )  =  si(n)  4-  S2(n).  Then 


(75)  = 

where  Wsii32(n ,  6)  is  the  cross- Wigner  distribution  between  $i(n)  and  S2 (n).  In 
time-frequency  parlance,  the  terms  W31  (n,  6 )  and  W52(n,  6)  are  called  auto-terms  and 
the  term  25ft  ( W81  yS2  (n,  0))  is  called  a  cross- term. 

Suppose  s\  (n)  and  S2(n)  are  two  discrete-time  chirp  signals,  and  consider  the  type-II 
Wigner  distribution  of  s(n)  —  5i(n)  +  S2(n).  The  large  maxima  of  W*1  ,  9)  with 

respect  to  6  will  not  all  correspond  to  auto-terms;  some  of  them  will  correspond  to 
cross-terms.  However,  the  auto-terms  may  be  distinguished  from  the  cross-terms  as 
described  below. 

In  the  time-frequency  plane,  the  auto-terms  are  localized  to  lines  and  the  cross-terms 
are  located  half-way  between  the  auto-terms.  Although  the  cross-terms  can  be  as  large 
as  the  auto-terms,  they  exhibit  an  oscillatory  behaviour  in  the  direction  of  the  line  that 

4When  we  do  this,  the  meanings  of  Hq  and  Hi  will  change,  but  this  is  only  of  pedagogic  significance. 
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is  half-way  between  the  lines  to  which  the  auto-terms  are  localized.5  The  auto-terms  do 
not  exhibit  this  oscillatory  behaviour.  Therefore,  the  auto-terms  can  be  distinguished 
from  the  cross-terms  by  performing  line-integrations  in  the  time-frequency  plane. 
Equivalently,  we  can  perform  correlations  as  described  below. 

Whether  a  maximum  of  Wj1  0),  at  say  9  —  Cm,  corresponds  to  an  auto-term  or  a 
cross-term  can  be  decided  by  the  rule 


(76)  max  I^c^ca)  J  7 o, 

for  some  appropriate  70,  where  H auto  and  H cross  are  the  ‘auto-term’  and  ‘cross-term’ 
hypotheses.  If  an  auto-term  is  detected  then  the  above  maximization  procedure  will 
yield  its  frequency-rate  parameter. 

Therefore,  in  the  presence  of  noise,  multiple  chirp  signals  can  be  detected  by  using  the 
following  two-stage  procedure: 

1 .  Identify  all  local  maxima  of  W*1  ,  9)  with  respect  to  9  that  exceed  7m, 

2.  For  each  local  maxima  identified  in  step  1,  perform  this  test:  Suppose  0  =  Cm  is  an 
identified  local  maximum,  then  maximize  rr(Cm,C2)  with  respect  to  02.  If  the 
maximum  exceeds  /ya  declare  the  presence  of  a  chirp  signal  with  median  frequency 
Cm  and  frequency -rate  c<i. 

Note  that  this  two-stage  procedure  is  not  different  from  the  one  described  for  detecting 
single  chirp  signals. 


5.3  Further  Enhancements 

When  considered  as  a  function  of  0,  WsI(^f^y0)  will  have  global  maximum 
(main-lobe)  at  the  median  frequency  and  local  maxima  (side-lobes)  on  both  sides.  The 
ratio  between  the  global  maximum  and  the  next  largest  local  maximum,  known  as  the 
peak  sidelobe  ratio,  is  about  13  dB.  Thus  a  target  can  get  masked  by  another  target  in  its 
vicinity  that  is  larger  by  13  dB  or  more.  To  overcome  this,  we  need  to  use  some 
modified  version  of  W/J(^^,0)  obtained  by  using  a  weighting  function  before  taking 
the  Fourier  series  sum. 

The  location  of  the  cross-terms  have  a  logical  relation  to  the  location  of  the  auto-terms. 
A  sophisticated  method  of  distinguishing  the  cross-terms  from  the  auto-terms  that 
makes  use  of  this  information  needs  to  be  derived. 

5Thc  exact  structure  of  a  cross-term  will  be  described  later. 
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Time-frequency  distributions  that  have  the  same  periodicity  properties  as  the  type-II 
Wigner  distribution  but  have  reduced  cross-terms  need  to  be  derived. 

Having  detected  the  real  chirp  signals  and  estimated  their  parameters,  we  can  refine 
those  estimates  by  using  them  as  starting  points  for  an  iterative  method  of  maximizing 
Ar(ci,C2).  So,  an  iterative  method  of  maximizing  Ar(ci,  C2)  needs  to  be  derived. 


5.4  Cost  of  Computations  Involved 

Assuming  we  have  some  prior  knowledge  of  the  range  of  frequency-rates  likely  to  be 
encountered,  the  cost  is  dominated  by  the  cost  of  computing  Wj1 

Recall  from  Section  4  that  for  N  even 


/  AT  —  2  \  1 

(77)  W*1  l  ,  6  J  =  Y,  r(N/2  +  k)r*(N/ 2  -  1  -  k)e~^k+^6. 

'  '  k——Nj2 

Obviously,  this  can  be  computed  at  M  uniformly  spaced  values  of  0  over  [— 7r/2,  7t/2] 
in  a  straightforward  way  by  using  the  Fast  Fourier  Transform  algorithm  provided 
M  >  N  and  M  is  an  integer  power  of  2.  Denoting 

a(k )  =  r(N/2  +  k)r*(N/ 2  —  1  —  k),  we  see  that  a(k)  —  a*(—(k  +  1)).  This 
symmetry  property  allows  us  to  reduce  the  cost  roughly  by  a  factor  of  two. 

5.5  Consequences  for  Tracking 

As  was  shown  in  Section  2,  accelerating  targets  give  rise  to  discrete-time  chirp  signals. 
Thus  we  can  use  the  method  described  in  this  section  to  detect  a  constantly  accelerating 
target  and  estimate  its  median  velocity  and  acceleration. 

As  far  as  tracking  a  target  is  concerned,  the  availability  of  some  form  of  velocity 
measurement  should  enhance  the  tracking  performance  over  tracking  based  on  range 
and  azimuth  measurements  alone.  A  tracking  algorithm  for  accelerating  targets  based 
on  range,  azimuth,  and  median  radial  velocity  measurements  should  be  easy  to  derive, 
especially  if  we  have  derived  a  tracking  algorithm  for  non-accelerating  targets  based  on 
range,  azimuth,  and  radial  velocity  measurements.  In  the  same  vein,  an  implementation 
of  a  tracker  for  non-accelerating  targets  based  on  range,  azimuth,  and  radial  velocity 
measurements  should  be  easily  upgradable  for  tracking  accelerating  targets  based  on 
range,  azimuth,  and  median  radial  velocity  measurements. 

To  make  use  of  the  acceleration  measurement,  in  addition  to  the  median  velocity 
measurement,  an  entirely  new  tracking  algorithm  needs  to  be  derived. 
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6.  Conclusion 


A  variety  of  methods  have  been  proposed  for  the  detection  of  a  signal,  with  unknown 
signal  parameters,  in  a  noisy  envirnoment.  In  this  report,  the  noise  statistics  are 
incorporated  to  reveal  that  certain  processing  of  the  Wigner  distribution  signal 
representation  can  lead  to  an  optimal,  and  often  easy  to  compute,  detection  scheme.  For 
the  special  case  of  linear  FM  signals  in  complex  white  Gaussian  noise,  it  is  shown  that 
the  optimal  detector  is  equivalent  to  integrating  the  Wigner  distribution  along  the  line  of 
instantaneous  frequency.  If  the  position  and  sweep  rate  of  the  linear  chirp  are  unknown, 
then  a  Generalized  Likelihood  Ratio  Test  (GLRT)  leads  one  to  integrate  the  Wigner 
distribution  along  all  possible  lines  in  the  time-frequency  plane  and  choose  the  largest 
integrated  value  for  comparison  to  a  threshold.  Simulation  examples  of  the  Wigner 
distribution  scheme  are  given  to  demonstrate  the  utility  concerning  the  detection  of  the 
proposed  method. 

In  this  report,  a  new  discrete-time  Wigner  distribution  has  been  also  proposed.  The 
ideas  presented  can  be  extented  to  formulate  a  Wigner  distribution  in  both  time  and 
frequency,  for  discrete-time  signals  of  finite  duration,  in  the  same  spirit  as  development 
of  discrete  Fourier  transform.  Application  of  these  ideas  to  actual  data  analysis  is  the 
subject  of  further  research. 
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Annex  A 

Proof  of  Equation  5.3  of  Section  2 


Recall  from  Section  3  the  definition 


(A.1) 


Ar(ci,ca)  = 


N—l 


Y  r(n)e~j(-cin+2C2n2) 


n— 0 


We  can  write  this  as  the  double  summation 


(A.2) 

(A.3) 


Ar(ci,C2) 


N—l  N—l 

EE  r(ni)r*(n2)e-^cl(m-n2>f2C2(n?-ni)>, 

Til— 0  722=0 


N- 1  N—l 

Y  Y  r(n1)r*(n2)e-j(ni-n2)ici+*C2(-ni+n2)). 

ni— 0  «2=0 


Next  we  break  the  double  summation  into  two  double  summations  -  one  double 
summation  being  over  m  and  n 2  where  (m  —  112)  is  even  and  the  other  double 
summation  being  over  m  and  n: 2  where  (ni  —  n2)  is  odd.  Thus  we  define 


(A.4)  A '(<*,<*)=  El'iEto  r,(ni)r*(n2)e_^ni_n2^cl+2C2(ni+n2^, 
(ni  —  ra2)  is  even 


and 


(A.5)  A"(c1,c2)=  En“oES  r(m)r*(n2)e-^-^+^<”1+112)). 
(ni  —  n2)  is  odd 

To  compute  A*  (ci,  c2),  we  use  the  variable  change  6 


(A.6)  ni  +  n2  =  2m, 

(A.7)  ni  —  n2  =  2k, 

6Note  that  n,i  +  m  is  even  if  and  only  if  ni  —  ni  is  even  for  integers  m  ,ni. 
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and  obtain 


TV— 1  min(m,AT— 1— m) 

(A. 8)  AIr(c1,c2)=  Y  r(m  +  k)r*(m-k)e-j2k(cl+C2Tn\ 

m—0  k—~  min(m,AT— 1—  m) 

TV— 1 

(A.9)  =  y;  WrJ(m,  Cl  +  c2m) 

m—0 
AT— 1 

(A- 10)  =  y  W?(m,0)|, 

=ci+C2m  * 

m=0 


Similarly,  to  compute  A^(ci,  C2),  we  use  the  variable  change 


(A.l  1)  Tt\  ~b  ^2  —  2tti  +  1, 

(A.12)  ni—ri2  —  2k  +  1, 


and  obtain 


N- 2  min(m,AT— 2— m) 

A"(d,ca)  =  Y  Y  r(m  +  /c  +  l)r*(m  —  fc)e-i(2fc+1)(cl+2C2+C2m)? 

m—0  /c——  min(m+l,TV- 1— m) 

(A.  13) 

AT-2 

(A.  1 4)  =  y  W?1  (m,  Cl  +  -  c2  +  c2m) 

m—0 

TV-2 

(A,i5)  = 

m=0 

By  combining  these,  we  obtain 


4A(t#)c2)=  A^(ci,c2)  +  A^(ci,c2), 

AT-1  TV-2  1 

(A. 17)  =  y  W/(m,  ci  +  c2m)  +  y  Wr/;  (m,  ci  +  -c2  +  c2m), 

m—0  m—0 

TV— 1  1  AT— 2 

=  y  W///(2m, ci  +  -c22m)  +  y  Wj.11  (2m  +  1,  Ci  +  -c2(2to  +  1)), 

m=0  m—0 


38 


DREO  TM  2001-141 


(A.  18) 


2N-2  1 

=  ^  W*II(m,c1  +  -c2m)  + 

m  =  0 

(m  is  even) 

2N-2  1 

(A.  19)  £  ^Kd  +  ^cam), 

m  =  0 
(m  is  odd) 

2N—2 

(A.20)  =  ^  +  -c2m), 

m=0 

2N-2 

(A.21)  =  2 

m— 0 

Thus  we  have  proved  equation  (28)  of  section  4.2. 
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Annex  B 

A  Historical  Review  of  the  Use  of  Wigner 
Distribution  for  Signal  Detection _ 

In  this  appendix,  we  provide  a  short  historical  review  of  the  use  of  the  Wigner 
distribution  for  signal  detection  in  additive  white  Gaussian  noise.  This  may  explain 
how  the  error  pointed  out  in  Section  4.4. 1  may  have  been  committed. 

B.1  The  Case  of  Continuous-Time  Signals 

For  a  conti nous-time  signal  r(t),  the  Wigner  distribution  is  defined  as  [13] 

(B.l)  WT(t,u))  =  J  r(t  +  r/2)r*(t  —  T/2)e~^Tdr. 

Suppose  we  have  observed  a  continuous- time  signal  r(t)  and  want  to  detect  the 
presence  or  absence  in  r,(f)  of  a  deterministic  signal  s(t)  with  the  backgorund  being 
additive  white  Gaussian  noise.  The  classical  maximum-likelihood  method  can  be  based 
on  the  rule 

Hx 

(B.2)  |/ r{t)s*{t)dt 

Ho 

where 

•  H0  is  the  Noise-Only  Hypothesis  -  that  is  r(t.)  —  w(t), 

•  Hi  is  the  Signal-plus-Noise  Hypothesis  -  that  is  r(t)  =  s(t )  +  w(t), 
where  w(t )  is  white  Gaussian  noise.  For  the  chirp  signal 

(B.3)  s(t)  =  e}'(W0t+2mt2), 

the  rule  is 

Hj 

(B.4)  j/ 

H0 
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In  [12]  the  above  rule  was  shown  to  be  approximately  equivalent  to 


a 


Hi 

/oo  > 

Wr(t,u>o+mt)  7 

OO  ^ 

H0 

for  chirp  signals  of  large  duration.  The  rule  was  rigourously  shown  to  be  valid  even  for 
chirp  signals  of  finite  duration  in  [16].  More  Specifically  it  was  shown  in  [16]  that 


2 

(B.6)  J  r(i)e-i(wot+amta>dt  =  J  Wr{t,u0 +mt)dt. 


A  more  comprehensive  study  of  the  use  of  the  time-frequency  distributions  for 
detecting  signals  is  given  by  [17]. 


B.2  The  Case  of  Discrete-Time  Signals 

For  a  discrete-time  signal  r(n),  the  Wigner  distribution  defined  by  [5]  has  become 
widely  accepted.  Their  definition  of  the  Wigner  distribution  is 


(B.7)  W?M(n,0)  =  2  J]V(n  +  k)r*(n  -  k)e~j2k0, 

k 

which  is  twice  the  type-I  Wigner  distribution  which  we  have  provided  in  Section  4.1.1. 

Suppose  we  have  observed  a  discrete-time  signal  r(n)  and  want  to  detect  the  presence 
or  absence  in  r(n)  of  a  deterministic  signal  s(n)  with  the  background  being  additive 
white  Gaussian  noise.  The  classical  maximum-likelihood  method  can  be  based  on  the 
rule 


% 


(B.8) 


2 


J~V(n)s*(n) 


Hx 

> 

< 

Ho 


where 

•  H0  is  the  Noise-Only  Hypothesis  -  that  is  r(n )  =  v(n), 

•  Hi  is  the  Signal-plus-Noise  Hypothesis  -  that  is  r(n )  =  s(n)  +  v(n), 
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where  v(n)  is  white  Gaussian  noise.  For  the  chirp  signal 


(B.9)  s(n)  =  eJ'(cin+2C2t,2>, 

the  rule  is 


(B.10) 


2 


^r(n)e-^Cin+2C2n2> 


n 


Hj 

> 

< 

Ho 


The  equivalence  of  equation  (104)  for  continuous-time  signals  was  erroneously 
assumed  to  carry  over  to  discrete-time  signals  as 


(B.  11) 


2 


^r(n)e"^cin+2C2n2) 


+c2n), 

n 


and  it  was  claimed  that  maximum-likelihood  detection  can  be  based  on  the  rule 

Hj 

(B.  12)  ^WfM(n,Cl-l-c2n)  >  7. 
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Annex  C 

The  Wigner  Distributions  of  a  Discrete-Time  Chirp 
Signal _ 


In  this  appendix  we  derive  the  type-I,  type-II,  and  type- HI  Wigner  distributions  for  the 
discrete-time  chirp  signal 


(C.l) 


J  gri&m+jfen2) 

l  o 


if  0  <  n  <  (N  —  1) 
otherwise. 


C.l  Type-1  Wigner  Distribution 

Recall  that  the  type-I  Wigner  distribution  W/(n,  0)  is  defined  as 


(C.l)  Wl (n,  0)  =  *(»  +  k)s*(n  ~  k)e~j2k9, 

k 

where  n  is  integer- valued  and  0  is  real- valued. 

It  is  easy  to  see  that  the  signal  product  term  s(n  +  k)s*(n  —  k )  is  zero  outside  the 
ranges  0  <  n  <  (N  —  1)  and  —  min(n,  N  —  1  —  n)  <  k  <  min(n,  N  —  1  —  n).  Thus 
Wl (n,  0)  =  0  outside  0  <  n  <  (N  —  1). 

Denote 


{n  if  n<  (N  —  l)/2, 

(N  —  l)/2  if  iV  is  odd  and  n  —  (N  —  l)/2, 
N-l-n  ifn>(N  —  l)/2. 


In  terms  of  this,  for  0  <  n  <  (N  —  1), 


* 


» 


(C.4) 

(C.5) 

(C.6) 


i 

Wj(n,0)  =  ^  ej(bl(n+k)+±b2(n+k)2)e-j(b1(n-k)+ib2(n-k)2)e-j2ke, 

k~—l 

l 

_  ej(6i2fc+^M«fc)e-j2A^^ 

k=-l 

l 

—  ^  ^  q— j(0—  (bi-\-b27i))2k 

k=-l 
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which  is  a  sum  of  a  geometric  series  that  can  be  easily  evaluated. 


We  first  observe  that  the  sum  can  be  written  as 


i 

(C.7)  e~jak * 

k=-l 

where  a  =  2(0  -  (&i  +  hn))  and  then  evaluate  it  to  be 


i 


(C.8) 


Y  e~jQk  = 


k=—l 

Thus,  for  0  <  n  <  (N  —  1), 


{ 


2/  +  1 
8in[o(l+§)] 
sm(f ) 


if  a  =  0  mod  27r, 
otherwise. 


j  (  2  min(n,  N  —  1  —  n)  +  1  if^  =  6i+  b^n  mod  n, 

(C.9W3  (n,0)  =  <  3in[(g-(fci+&2n))(2  min(n,AT-l-n)+l)l  otherwiSe 

^  sin[0-(6i+62n)] 


C.2  Type-ll  Wigner  Distribution 

Recall  that  the  type-II  Wigner  distribution  0)  is  defined  as 


(C.10)  W”(n,6)  =Ys(n  +  k  +  1)s*(^n  ~  k)e~i(2k+1)9, 

k 

where  n  is  integer- valued  and  0  is  real- valued. 

It  is  easy  to  see  that  the  signal  product  term  s(n  +  fc  +  l)s*(n  —  k)  zero  outside  the 
ranges  0  <  n  <  (N  —  2)  and  —  min(n  +  1,  N  —  1  —  n)  <  k  <  min(n,  N  —  2  —  n). 
Thus  W^iriyO)  =  0  outside  0  <  n  <  (N  —  2). 

Denote 


(C.  1 J  )=  min(n,  N  —  2  —  n) 


n  if  n  <  (N  —  2)/2 , 

i  (N  —  2) /2  if  N  is  even  and  n  =  (N  —  2) /2, 
k  N-2-n  if n>  (N  —  2)/2. 


In  terms  of  this,  for  0  <  n  <  (N  —  2), 


44 


DREO  TM  2001-141 


W^I(n,  6)  =  e?(fci(n+Ar+l)+^62(n+fc+l)2)e-j(fe1(n-fe)+^62(n-/c)2)e-j(2Aj+l)0^ 

fe=-(/+l) 

(C.12) 

i 

(C  13)  —  ^  ^  £?(&i  (2^+1)+ 2  ^2 (4nk+2n+2k+l))g--j(2k+l)$ 

fc=— (i+1) 

/ 

—  ^  ^  g— j(0— (fcl+|62+&2^))(2A:+l) 

fc— —  (i+1) 


(C.  14) 

(C.15) 


=  e-j(^-(ti+|62+62n)) 


e~j  (0-  (&i +|  fo-{-b2n))2k  ^ 

fe=-(i+ 1) 


which  is  a  scaled  version  of  a  sum  of  a  geometric  series  that  can  be  easily  evaluated. 
We  first  observe  that  the  sum  can  be  written  as 

i 

(C.16)  ]T  e~jak, 

k=-(l+ 1) 


where  a  =  2(9  —  (&i  +  |&2  +  &2^))  and  then  evaluate  it  to  be 


(C.  17) 


l 


Y  e~jak 

k=-(i+i) 


if  a  =  0  mod  27T, 


Now,  introducing  the  scaling  factor,  we  see  that 


(C.  18) 

(C.19) 


i 


k=—(l+ 1) 


e  *  %  (21 +  2) 

sin  [q(  1+1)] 
sin(f)  • 


if  a  =  0  mod  27t, 


(21  +  2) 

-(2Z  +  2) 

sin[q(f+l)1 
sin(f)  ' 


if  |  =  0  mod  2n, 
if  |  =  7 r  mod  27r, 


Thus,  for  0  <n<(N  —  2), 
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w/'M) 

(C.20) 


2  min  (n,  N  —  2  —  n)  +  2 
—(2  min(n,  AT  —  2  —  n)  +  2) 

sin[(fl—  (6i  +  §62+^2^))(2min(nrAT— 2-n)+2)] 
sin  [0 — (&i  +  2  ^2 + ^2n)] 


if  0  —  ?>i  +  t|&2  +  mod  27r, 

if  0  =  7T  +  fei  +  ^&2  +  mod  27T, 

otherwise. 


It  is  useful  to  write  this  as 


w/'M) 

(C.21) 


2mm(n  + 


t  2 > — j.  —  \,v*2  ))' 

2min(n  +  ^,N  -  1  -  (n+  ^))  +  ll 

(0-  (bi  +62  (n+ 1 )))  (2  min(n+ ^ ,  N-  l-(n+ \ 

.  I*  /i  .  .  /  .1 


»+0 


sin[g- (61 +62  (n+ 1 ))] 


if  8  =  bi  +  b2(n  +  |)  mod  2n, 
if  0  =  7r  +  bi  +  62 (n  +  |)  mod  27r, 

otherwise. 


C.3  Type-Ill  Wigner  Distribution 

Recall  that  the  type-III  Wigner  distribution  Wf11  (n,  0)  is  defined  in  terms  of  the  type-I 
and  type-II  Wigner  distributions  as 


(C.22) 


W/(»/ 2,0) 
Wl\{n-  \)M 


for  even  n, 
for  odd  n. 


For  n  even  and  0  <  n  <  2 (N  —  1), 


(  2min(^, N  —  1  —  |)  +  1 
(<W#j(n,0)  =  i  sinf(0-(6,  +  i62n))(2min(f,tV-l-f)+l)l 

I  - sin[9_(Va  +  [fr2n)] 


if  0  =  61  +  \b2r1  mod  7T, 
otherwise. 


For  n  odd  and  0  <  n  <  2 (N  —  1), 


w/"M) 

(C.24) 


2min(f,JV-l-f)  +  l 
-  2min(s,AT-l-f)  +  l] 

sin  (g—  (fci  +  ^ &2n))(2min(^,iy-l  —  § )+l)1 


sin[0-(6i+|62n)J 


if  8  =  bi  +  \  b^n  mod  2n, 
if  0  =  7r  +  61  +  \b2r1  mod  27r, 

otherwise. 


Combining  the  above,  for  all  0  <  n  <  2(iV  —  1), 
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Wain(n,6)  =  | 
(C.25) 


2mm(f,jV-l-f)  +  l 
(-1)"  [2min(f,iV-l-f)  +  l] 

sin[(9-(fei+jfen))(2min(g,Ar-l-g)+l)] 

sin[@-(6i+|t>2n)] 


if  0  =  b\  + fan  mod  2  it, 
if  6  =  7r  +  &i  +  fan  mod  2n, 

otherwise. 


4 


« 


i 
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